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We study the critical behavior of the chiral condensate near the QCD phase transition in the background of 
two fixed light dynamical (sea) quarks. We study the condensate for 5.245 < /3 < 5.3 and 10"^" < mvai < 10 
(in lattice units) on a 16'^ x 4 lattice using staggered fermions with rrisoa fixed at 0.01. 



1. INTRODUCTION 

Simulations of the QCD chiral phase transition 
with dynamical quarks have been carried out in 
the past for a variety of quark masses at differ- 
ent lattice spacings[Q. For two light flavors of 
staggered fermions these simulations indicate a 
smooth transition. If this turns out to be a sec- 
ond order transition, one expects the associated 
critical behavior to be in the 0(4) universality 
class [Ij. There have been attempts [|| to fit the 
existing data with the 0(4) exponents with some 
success. However to gain confidence in such fits 
one must understand the physics better, espe- 
cially in the range of parameters used in these 
simulations. For example, when the quark masses 
are sufficiently large, a weak first order signal may 
not be resolved and could be seen as a smooth 
transition. In order to answer this question and 
understand the physics of the transition from a 
slightly different angle we have asked a somewhat 
new question; how does the transition look when 
studied in the small mvai (valence-quark mass) 
limit, for a small but fixed mgca (sea-quark mass)? 

2. CHIRAL CONDENSATE 

The order parameter of QCD is the chiral con- 
densate (i/'V') which is given by the formula 
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Here D^u] is the Dirac operator on the lattice 
for the gauge configuration [U] and m the quark 
mass in lattice units. The expectation value is 
taken over the distribution of gauge configura- 
tions given by the usual QCD path integral 



(O) 



d[U]0[U] exp[-/3Sg + Sf(m)] (2) 



where Sg is the Wilson gauge action with /3 = 
6/(7^ and exp[Sf (m)] is the staggered fermion de- 
terminant for two quark flavors of mass m. In our 
convention, is normalized to \/m for large 

m by choosing fi. 

2.1. Partial quenching 

Distinguishing the masses that enter (1) im- 
plicitly through the fermion determinant and ex- 
plicitly through the trace of the propagator, one 
gets 
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In this study we compute (■0V'('^vai))/3,msca fo^' 
various values of mvai at a small but fixed value 
of TTisoa near the critical coupling /3c. This is rel- 
atively easy in comparison to changing the val- 
ues of mvai and mgea simultaneously as required 
in full QCD. However such studies are familiar 
at zero temperature (far below the critical cou- 
pling) and is referred to as the partially-quenched 
approximation. In this work we have extended 
such a study to finite temperatures especially to 
get a better understanding of the phase transi- 
tion. ('!/'^("ivai))/3,ms„a can be referred to as the 
partially-quenched chiral condensate. 
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2.2. Dirac spectral density 

If one defines P[£/](A) as the density(per unit 
volume) of eigenvalues of the anti-hermitian Dirac 
operator D^uj at the eigenvalue iA, then one can 
rewrite (3) as 



('0-0(mval))/3,ms. 
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Wsea — * limit we obtain the 
Casher formula H 
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When p^o(O) 7^ we have chiral symmetry 
breaking. Furthermore p^,„^^^(0) is expected to 
be a regular function of irises for small rrisca at 
least far from the critical coupling where chiral 
symmetry restoration takes place. Hence we ex- 
pect p^^m^^^(O) 7^ below the critical coupling 
and zero above. This suggests that m^^S'^^ ^^^^ 
show a striking critical behavior around the crit- 
ical coupling f3c. This critical behavior can be 
observed in (V'V'(™vai))/3,msoa due to (4). 

3. CRITICAL BEHAVIOR 

In this study we have computed 
(-0'0('7ivai))/3,ms<=a using (3) at a fixed msea = 0.01 
in a spatial volume of 16'^ and Nt = 4. We have 
studied (V'V'(»7tvai))/3,m,,a for 10~^° < mvai < 10 
which allows us to easily identify the finite vol- 
ume and finite lattice spacing regions. We have 
simulated at nine different couplings in the region 
5.245 < P < 5.3. Figure ^ shows some of these 
results on a log-log plot. 

3.1. Finite lattice effects 

The range of Wvai used in our study gives us an 
idea of the finite lattice-spacing and finite volume 
effects. We see that for mvai > 0.01 finite lattice 
spacing effects become important and for m^ai < 
2.5 X 10~^ finite volume effects begin to set it. 
One knows that in a finite volume for sufficiently 
small masses the condensate must go linearly to 
zero. This is convincingly seen in Figure |^ for 
sufficiently small masses. The region 2.5 x 10^'* < 
Wvai < 10^2 could be taken to be free of the 
finite size and finite lattice spacing effects. This 
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Figure 1. A log-log plot of (V'V') vs. mvai for dif- 
ferent values of /3 at a fixed mgca = 0.01. The 
figure clearly shows finite volume and finite lat- 
tice spacing effects. 



Sf:S^.ol]pi§^*^^^t-!&i%efl5f of couplings around 
aeT=ffiiMtion one expects 
a pure power law of the form {^pip) ~ ^A, the 
critical point. It is interesting to see if such a form 
also describes the dependence of (V'V'('7tvai))/3,maea 
as rrivai — > 0. With the available statistics we 
find that using a simple least squares fit the data 
is consistent with a power law over a range of 
couplings with the power j ranging from about 
0.02 at a /3 of 5.245 to 0.8 at a /3 of 5.3. JThese fits 
are shown on Figure I The plot of {i^iP) vs. f3 
for a range of mvai is given in Figure ^. Here we 
again sec evidence for a smooth transition with 
ftle?MI^fil^5i^ffi|^tween 5.26 and 5.275. 

The data suggests that at mgoa = 0.01 there is 
no discontinuity in the distribution of gauge con- 
figurations as a function the coupling. Any possi- 
ble discontinuity will be enhanced by our method 
of measuring the condensate at very small valence 
quark masses. Hence this study supports our ear- 
lier results that at msea = 0.01 there is no evi- 
dence of a first order transition. Given the accu- 
rate results that we have obtained for msea ~ 0.01 
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close to the critical region we can compare with 
our earlier results for msca = 0.025. Figure ^ is a 
plot of the entropy density from the present run 
and the earlier run at mgoa of 0.025. Except for 
a change in Pc caused by the variation of Wsea, 
the sharpness in the two graphs appear remark- 
ably similar unlike the graphs for m-vai = 0.01 and 
0.025 in Figure H where msea = 0.01 is fixed. This 
comparison of Figure ^ and Figure ^ also shows 
that the sea-quark effects are the dominant effects 
in the analysis of pseudo-critical coupling done 
by Karsch|^. Studying the effects of msca and 
nival separately will be helpful in understanding 
some of these observations and can make partially 
quenched studies a useful technique in learning 
about the real critical region. 
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2. The plot of l^t/^} vs. mvai is con- 
with power-law behavior over a range of 
gs, in the region free of finite lattice ef- 
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Figure 3. The plot of {ip'ip) ^s. (3 for various val- 
ues mvai at a fixed mgca = 0.01 shows a smooth 
transition. 
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Figure 4. The plot of entropy density vs. (3 at 

rrisea = "^vai of 0.01 and 0.025. 



